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We study water adsorption-induced deformation of a monolithic, mesoporous silicon membrane tra¬ 
versed by independent channels of ~8 nm diameter. We focus on the elastic constant associated with 
the Laplace pressure-induced deformation of the membrane upon capillary condensation, i.e. the 
pore-load modulus. We perform finite-element method (FEM) simulations of the adsorption-induced 
deformation of hexagonal and square lattices of cylindrical pores representing the membrane. We 
find that the pore-load modulus weakly depends on the geometrical arrangement of pores, and can 
be expressed as a function of porosity. We propose an analytical model which relates the pore-load 
modulus to the porosity and to the elastic properties of bulk silicon (Young’s modulus and Poisson’s 
ratio), and provides an excellent agreement with FEM results. We find good agreement between our 
experimental data and the predictions of the analytical model, with the Young’s modulus of the pore 
walls slightly lower than the bulk value. This model is applicable to a large class of materials with 
morphologies similar to mesoporous silicon. Moreover, our findings suggest that liquid condensation 
experiments allow one to elegantly access the elastic constants of a mesoporous medium. 


Mesoporous silicon (pSi) prepared by electrochemical 
etching of bulk silicon has been attracting much atten¬ 
tion from both fundamental and applied sciences owing 
to its unique optical, electrical and thermal properties 
m- While control over mechanical properties of pSi 
is necessary for its applications, the understanding of its 
response to mechanical loads is mainly limited to mea¬ 
suring Young’s modulus of porous samples dlHl. 

Since pSi can be prepared as a bulk (monolithic) meso¬ 
porous system with parallel, channel-like independent 
pores, it has been of a particular interest for studies of 
fluid adsorption On the one hand, this matrix al¬ 
lows one to study the influence of spatial confinement 
on the physics and chemistry of liquids PfT5]. On the 
other hand, it is possible to synthesize composite materi¬ 
als with finely tuned optical and electrical properties by 
liquid adsorption or infiltration mile]. Therefore, explo¬ 
ration of the mechanical properties of silicon relevant to 
interactions with liquids, the topic of this Letter, is also 
of high importance. 

When a fluid is adsorbed in a mesopore, it exerts a 
pressure on the pore walls, which is typically of the order 
of 10^ Pa [T7]. This pressure causes deformation of the 
pore, and as a result deformation of the porous material 
as a whole. This effect, known as adsorption-induced 
deformation has been experimentally observed for var¬ 
ious mesoporous materials: Vycor glass na El] , tern- 
plated silica |2nU2^ . low-k films [SSIEl], aerogels [^l26] . 
porous gold m and pSi |S5|- There are two ways to 
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measure the adsorption strains experimentally: for ma¬ 
terials which can be prepared as macroscopic samples, 
the dilatometric technique can be used, and the reported 
strain is the relative change of the length of the sample 
[m iiiiesiet]. For crystalline materials or materials 
which have periodic pore structure. X-ray diffraction can 
be used, and the strain is calculated as a change of the 
crystal da HI] or pore-spacing lattice constant dniEi]. 
Irrespective of the technique, the measured strains for all 
these materials (except for aerogels) are of the order of 
10“'* — 10“^. Such small strains are in the linear elastic 
regime, and it is reasonable to assume a linear relation 
between the pressure in the pore and the experimentally 
observed strain with a proportionality constant M, called 
the pore-load modulus [H] . If the pore-load modulus for 
a certain porous material is known, it is possible to cal¬ 
culate the fluid pressure in the pores based on the ex¬ 
perimental data on adsorption-induced strain, providing 
information on the thermodynamics of the confined fluid. 
Also, understanding of elastic response of a porous ma¬ 
terial to adsorption is useful for its application to sensing 
and actuation 

The pore-load modulus relates the pressure inside the 
pore to the overall deformation of the sample, rather than 
relating an external load to a deformation, like Young’s or 
bulk moduli. Therefore an important question is how to 
relate the pore-load modulus to the material properties 
of the matrix and the pore geometry. In general, for a 
material with a wide pore size distribution (PSD) and 
arbitrary pore morphology and orientation, this question 
may be complicated. However, for systems such as pSi 
considered in this work, where the geometry is regular, 
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it can be resolved. 

Here we present a dilatometric study of the defor¬ 
mation of a macroscopic pSi membrane induced by ad¬ 
sorption of water vapor. From the experimental strain 
isotherm we calculate the pore-load modulus. We per¬ 
form FEM simulations of adsorption-induced deforma¬ 
tion of samples with hexagonal and square lattices of 
cylindrical pores. We find that the pore-load modulus 
weakly depends on the geometrical arrangement of pores, 
and can be expressed as a function of porosity. Therefore 
we relate the deformation of a porous sample to the de¬ 
formation of a single cylindrical tube from the pressure 
applied to its inner surface. Based on this model we de¬ 
rive an analytical expression for the pore-load modulus 
as a function of the porosity and the elastic properties of 
the non-porous material (Young’s modulus and Poisson’s 
ratio). The predictions of our analytical model are in ex¬ 
cellent agreement with FEM results. We also achieve 
good agreement between our experiment and our ana¬ 
lytical model for adsorption-induced deformation of pSi, 
suggesting that the Young’s modulus of silicon pore walls 
could be only slightly lower than the Young’s modulus of 
bulk silicon. 

Adsorption-induced deformation of pSi has previously 
been studied by Dolino et al. [28]. They obtained the 
strain normal to the plane of the sample from the shifts 
of X-ray diffraction peaks of the crystalline lattice of the 
silicon matrix, and calculated the normal-to-the-plane 
pore-load modulus from a comparison of the experimen¬ 
tal strain to the fluid pressure in the gas pressure region, 
where the pores are filled with capillary condensate and 
the fluid pressure in the pore has a simple form. They 
also proposed an analytical expression for the pore-load 
modulus as a function of material parameters based on 
Scherer’s model for porous glasses m and the theory 
of elasticity for cellular materials [32]. However, their 
model did not agree with the experimental data. Re¬ 
cently Grosman et al. reported a dilatometric study of 
both the in-plane and normal-to-the-plane deformation of 
pSi induced by n-heptane adsorption |33j . In this work 
the authors calculated the pore-load moduli for two sam¬ 
ples, and proposed their own mechanical model to relate 
these moduli to the Young’s modulus of silicon. By ap¬ 
plying their model to their experimental data the authors 
came to the conclusion that the Young’s modulus of sil¬ 
icon walls is lower than the Young’s modulus of bulk 
silicon by a factor of five. 

Another attempt to relate the pore-load modulus to 
material parameters was made by Prass et al. [H]. They 
studied adsorption strains in MCM-41 and SBA-15 sil¬ 
ica, which have morphologies similar to pSi, and pro¬ 
posed that the strain, derived from the change of the pore 
lattice parameter, behaves similarly to the hoop stress 
in thin-walled cylinders. Although that model provided 
good agreement with their experimental data, it was not 
consistent with their finite element method (FEM) cal¬ 
culations. Our model does not employ the thin-wall as¬ 
sumption and is fully supported with FEM calculations. 


Mesoporous silicon was prepared by electrochemical 
etching of highly boron-doped (100) silicon wafers (pro¬ 
ducer: SiMat, Landsberg, Germany; specific conductiv¬ 
ity: p=0.01-0.025 B cm). After an etching depth of 30/rm 
a high voltage was applied in order to remove the meso¬ 
porous part from the bulk silicon underneath. The re¬ 
sulting parallel, non-interconnected, tubular pores have 
a short-range hexagonal order. A scanning electron mi¬ 
croscope (SEM) image of the sample is shown in Fig. 



FIG. 1. SEM image of the pSi structure. Two-dimensonal fi¬ 
nite element representation of a porous material consisting of 
parallel cylindrical pores distributed on a hexagonal lattice. 
The mesh shown here is much coarser than the mesh used 
for the FEM calculations. Our FEM analysis shows that the 
deformation of the porous sample can be reduced to the de¬ 
formation of the cylindrical domains in the Re vicinity of each 
pore; the latter problem can be solved analytically. 

The sample was cut in a rectangular shape with 4 
mm X 7 mm and a thickness of 30 /rm. The sample’s 
long axis was parallel to the crystallographic (010) di¬ 
rection. As prepared the inner pore walls are Si-H ter¬ 
minated. A 30 minutes infiltration of a peroxide/water 
solution and subsequent rinsing with Millipore water ren¬ 
dered the inner pore space Si-OH terminated, and thus 
hydrophilic. A volumetric nitrogen adsorption isotherm 
measurement performed at 77 K and analysed within the 
NLDFT model [M] indicated a mean pore diameter of 
8 nm (width of the PSD 15%) and a volume porosity of 
60%. 



FIG. 2. (color online) Schematic of the experimental setup 
used for measuring the adsorption-induced deformation of 
porous silicon. The sample holder (dark grey, on the left) 
is connected to a load cell, and the sample holder on the right 
is connected to a linear motor stage to measure the length 
change AL Thick arrows (yellow online) indicate the force on 
the pore walls due to Laplace pressure upon water condensa¬ 
tion. 

The measurements of the adsorption-induced strains 
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were performed using a custom-made dilatometric setup, 
Fig. [2] The samples were tested in a sealed chamber, kept 
at a constant temperature of 24° C by means of a water 
circulation thermostat (Huber). The humidity inside the 
chamber was controlled by means of a Wetsys (Setaram) 
humidity generator, which was working with a flow of 200 
ml/min. The samples were clamped to two aluminium 
holders. The macroscopic deformation was measured by 
one of the holders that connects to a Physics Instruments 
M-404 linear motor stage (resolution of 2 ^m), while the 
axial tensile force was measured using a Honeywell R- 
30 load cell (50 N max. load), attached to the other 
holder. The standard deviation of the measured force 
background noise over more than 100,000 points was 10 
mN. The strain of the sample was measured by driving 
the motors to keep a constant (small) force of 100 mN, 
while the humidity was continuously changed between 
5% and 95% RH at a rate of 10% RH/h. The stress 
on the sample under these conditions is estimated to be 
lower than 0.5 MPa. 

In Fig.|^we show the experimentally measured, macro¬ 
scopic strain of the sample a function of 

relative vapor pressure (= relative humidity) p/po = RH 
(po is the saturation pressure of water at 24°C). Here 
Al is the measured length change of the membrane and 
Iq the length of the sample at RH=0%. The behavior 
observed is typical of adsorption-induced deformations 
of mesoporous materials and is due to the change of the 
pressure P of the confined fluid (adsorption stress) in the 
mesopores The shape of the strain at low pressures 
(before the hysteresis loop related to capillary conden¬ 
sation) may differ and depends on the strength of the 
solid-fluid interactions [35]. At higher pressures, when 
the pores are filled with capillary condensate, the strain 
changes linearly with ln(p/po)- This trend is not specific 
to pSi and has been observed for all mesoporous materi¬ 
als investigated so far [I8ll26l |28l [33]. In our case, this 
region corresponds to pressures p/po ~ 0.83 — 0.93. The 
strain for this region is shown in the inset of Fig. 

When a pore is filled with fluid, the pressure P that 
the fluid exerts on a pore wall is given by HZ]: 


P = -^si/R + Pl, (1) 

where jsi is the solid-liquid surface energy, R is the pore 
radius and Pl is the Laplace pressure at the liquid menis¬ 
cus. The latter can be written as (Kelvin-Laplace equa¬ 
tion) 


Pl = 




Vi 



( 2 ) 


where Rg is the gas constant, T is the temperature and V) 
is the molar volume of the fluid. Since both the pressure 
in the pore and the experimentally measured strain of the 
porous sample as a whole vary linearly with ln(p/po)) it is 
easy to get the pore-load modulus from the following rela¬ 
tion P = Mem + C, where C is a constant, related to the 
surface energy of the solid-fluid interface. Although the 



FIG. 3. Measured strain of a monolithic mesoporous silicon 
membrane as a function of relative water vapor pressure p/po- 
Inset: Measured strain in the capillary condensation regime, 
where the strain depends logarithmically on p/po and thus 
linearly on the Laplace pressures in the pores. The solid line 
represents a linear fit of the strain as a function of log relative 
pressure, and its slope corresponds to the pore-load modulus 
M. 


value of C is important for studying the thermodynamics 
of confined fluids [35], it does not affect the present dis¬ 
cussion. Fitting of the data shown in the inset of Fig. 
gives M = 34.5 GPa. 

The key problem is to find the relation between the 
pore-load modulus M and the elastic constants of the 
non-porous material: Young’s modulus E and Poisson’s 
ratio V. We assume isotropic elasticity. Note that by the 
Young’s modulus E we mean the modulus of the pore 
walls, and not the effective Young’s modulus of a pSi 
sample as a whole {Ep), often reported in the pSi liter¬ 
ature [37]. These two moduli should not be confused, 
since the latter is a property of a solid-void composite 
and is a strong function of the porosity |7||38]. 

The deformation of a system with multiple pores can 
be solved numerically using FEM [33] . Since for the con¬ 
sidered material the PSD is narrow, the pores are parallel 
and distributed with a short-range hexagonal order, the 
pore-load modulus can be approximately calculated from 
deformation of a two-dimensional hexagonal lattice. The 
FEM representation of the problem is shown in Fig. 
The left and bottom surface are constrained in the x and 
y directions respectively, while the right and top surface 
move freely. The pressure P is applied in every pore, 
and the pore load modulus is estimated as the ratio of 
the applied pressure to the average engineering strain in 
the vertical direction. Poisson ratio v = 0.28 was used, 
corresponding to bulk silicon properties [30] • Systems 
of two different sizes, 5x8 (68 pores) and 10 x 16 (295 
pores), were used to check for finite size effects. We found 
that the results for the larger system differs from the re¬ 
sults for the smaller one by 2.5% at most. To evaluate 
the role of the lattice geometry we also performed the 
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simulations on a square lattice 14 x 14 (296 pores); the 
difference between the modulus calculated for square and 
hexagonal lattices is 4% at porosities (j) = 35% and 14% 
at (/> = 65%. 

The pore-load modulus calculated from FEM simula¬ 
tions is shown in Fig. We see that the moduli for 
different lattice geometries are close, and for a given E 
and v depend mainly on porosity (j) - the ratio of the 
volume of the voids to the total volume of the sample. 
For an analytical estimate of this dependence, we use the 
following considerations based on classical theory of elas¬ 
ticity. A single cylindrical pore with inner radius R and 
pressure P induces a purely deviatoric stress in an infinite 
plate, not leading to any swelling. Swelling is a conse¬ 
quence of the interaction of the pressurized pore with 
a free surface (where the normal stress must be zero). 
This is well known for the swelling of a pressurized cylin¬ 
der with outer radius Re (Fig. where the engineering 
strain is calculated from the increase of Re mma: 


_ u{Re) _ 

Re E{l-0 ’ 


(3) 


where ^ = R?/R^. In a large plate containing many 
pores arranged e.g., on a hexagonal or a square lattice, 
symmetry implies that the stress must vanish at half the 
distance between neighboring pores. Hence, the forces 
compensate at this point making it actually equivalent to 
a free surface. This suggests that the dilatational strain 
in the plate with many pores can be approximated by 
Eq. i where Re is taken as half the nearest neighbor 
distance in the pore lattice. The remaining question is 
how to relate ^ to (j). One possibility would be to use the 
lattice geometry to relate nearest neighbor distance to 
pore fraction, although this neglects the mechanical in¬ 
fluence of bits of material outside the symmetry planes. 
This would give the relation ^ = {2^/Z/'k)4> for the hexag¬ 
onal and ^ = (4/7r)(() for the square lattice, respectively. 
The best fit to the FEM data, however, is obtained with 
^ = 4> and, in the absence of a more detailed analyti¬ 
cal theory, we use this and write the pore-load modulus 
approximately as 

= 20 ^ ■ <“> 

Note that Eq. gives the the modulus in terms of the 
basic material properties E, v and porosity </). 

The results for the pore-load modulus of porous silicon 
are shown in Fig. The solid line shows the analyti¬ 
cal model Eq. which is in excellent agreement with 
the modulus M derived from our FEM calculations for a 
porous body. Fig. |^also shows that our model predicts 
a higher pore-load modulus than the model by Prass et 
al. [5T], especially at porosities below ~ 60%. Such de¬ 
viation is related to the thin wall approximation used in 
[21], which fails when the porosity decreases. 

Comparison of the model to our experimental data 
(Fig. 1^ shows that if we assume that the Young’s mod¬ 



FIG. 4. Ratio of the pore-load modulus to the Young’s mod¬ 
ulus of non-porous material as a function of porosity (j) (at 
V = 0.28). The results of FEM simulations of deformation of 
samples with hexagonal and square lattice show that the pore¬ 
load modulus weakly depends on the geometrical arrange¬ 
ment of pores, and can be expressed as a function of porosity. 
The proposed analytical model is in excellent agreement with 
FEM results. It also shows reasonably good agreement with 
the pore-load modulus calculated from the experimental data 
presented in this work. For comparison, the Young’s modulus 
of silicon was taken (100) E = 130 GPa | 40 |. 


ulus of pSi walls equals to the bulk (100) value, the ex¬ 
perimental data is close to the theoretical predictions, 
yet slightly below them. This deviation can be due to 
the irregularity of the pSi structure. Both the devia¬ 
tion of pore sizes from the mean size (i.e. PSD) and 
deviation of the pore arrangement from a perfect lattice 
lead to presence of thinner pore walls, which reduces the 
pore-load modulus of the structure. Also, the Young’s 
modulus of the walls can be somewhat lower than the 
modulus for bulk silicon; applying our model to calcu¬ 
late the Young’s modulus from our experimental value of 
M, we get E = 95 GPa. Note, however, that here the 
deviation of Young’s modulus of the pore walls from the 
Young’s modulus of bulk silicon is only 27% and not a 
factor of five as reported in [38j . 

The system of hexagonally-ordered parallel cylindrical 
pores is not unique to mesoporous silicon: a number of 
templated mesoporous materials have similar morphol¬ 
ogy, among them MCM-41 and SBA-15 silica, and silica 
monoliths with hierarchical pores, synthesized recently 
[22] . The model for the pore-load modulus proposed in 
this work can be applied to develop quantitative models 
of adsorption-induced deformation of these materials. It 
should provide better predictions than the earlier model 
m at porosities 60% and below, typical for mesoporous 
silica. 

To summarize, we have presented an experimental 
study of water adsorption-induced deformations of a 
monolithic, mesoporous silicon membrane, and extracted 
the pore-load modulus. We have proposed an analytical 
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model which relates the pore-load modulus to the poros¬ 
ity and to the elastic properties of bulk silicon, and jus¬ 
tified it by comparing with FEM simulations. We have 
found good agreement between our experimental data 
and the predictions of our model, with the Young’s mod¬ 
ulus of the pore walls slightly lower than the bulk value 
for silicon. Our model can be applied to a large class of 
materials with morphologies similar to mesoporous sili¬ 
con. Moreover, our findings suggest that liquid condensa¬ 
tion experiments allow one to access the elastic constants 
of a mesoporous medium. 
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